Astronomy & Astrophysics manuscript no. XUV5 REF BF5 


©ESO 2010 


June 25, 2010 





Research Note 

NGC 4262: a Virgo galaxy with an extended ultraviolet ring 

D. Bettoni 1 L.M. Buson 1 and G. Galletta 2 



o 

<N 

a 

<N 

o 

u 

6 

c3 



> 
r-- 

O 
O 



X 



1 INAF Osservatorio Astronomico di Padova, vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 
e-mail: daniela.bettoni@oapd.inaf. it , lucio.buson@oapd.inaf. it 

2 Dipartimento di Astronomia, vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 
e-mail: giuseppe . galletta@unipd . it 



Received ... ; accepted .... 2010 



ABSTRACT 



Context. The Galaxy Ultraviolet Explorer (GALEX) satellite has recently shown the presence of an extended, outer ring studded with 
UV-bright knots surrounding the lenticular galaxy NGC 4262. Such a structure — not detected in the optical — is coupled with a ring 
of atomic (HI) gas. 

Aims. We want to show that both star-forming and HI rings surrounding this SBO galaxy share the same radial distance from the 
galaxy center and spatial orientation. We want also to model the kinematics of the ring(s) and of the galaxy body. 
Methods. We make use of archive FUV and NUV GALEX data plus HI observations from the literature. 

Results. We confirm that the UV-bright and atomic gas rings of NGC 4262 have the same extent and projected spatial orientation. 
Their kinematics is not coupled with that of the galaxy stars. 

Conclusions. It is possible that NGC 4262 has undergone a major gas stripping event in the past which gave origin to the present 
"necklace" of UV-bright knots. 

Key words, galaxies: general — galaxies: elliptical and lenticular, cD — galaxies: individual: NGC 4262 — galaxies: interactions 



1. Introduction 

Before the advent of the Galaxy Evolution Explorer (GALEX) 
UV-bright galaxy features without counterpart in the optical 

■ were unavoidably neglected. In particular GALEX has shown 
that over 30% of spiral galaxies possess previously unknown 
UV-bright extensions of their optical disks (e.g. Thilker et al. 
2007). This is also the case of the lenticular galaxy NGC 4262 
in Virgo discussed here and surrounded by an extended, UV- 

■ emitting ring lying well beyond its optical image. 

As far as NGC 4262 is concerned, Hubble Space Telescope 
Advanced Camera for Surveys (HST/ACS) imaging (Cote et al. 
2006) unequivocally shows its morphology being typical of a 
SBO at optical wavelengths. More precisely, the galaxy has a 
bright, nearly spherical bulge, a bar, and a lens extending well 
outside the bar radius (Laurikainen et al. 2010). 

GALEX NUV (1771-2831 A) and FUV (1344-1786 A) im- 
ages (Fig. 1) show instead an external knotty feature, something 
like "beads on a string" likely made of individual, hot star clus- 
ters. In this respect, NGC 4262 — having clustered UV-bright 
sources in its outer parts — could be classified as a Type 1 ex- 
tended ultraviolet disk (XUV), a kind of feature discussed in de- 
tail by Thilker et al. (2007) and properly modeled by Bush et al. 
(2010); note, however, that the lack of H a observations does not 
allow us to state for certain we are looking at typical HII regions. 

In this paper we want to put the question whether the (strict) 
coincidence of the UV knots and HI ring does reflect a past in- 
teraction event. The origin of the UV-bright ring is of necessity 
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Fig. 1. GALEX combined FUV and NUV image of NGC 4262 
(blue=FUV, yellow=NUV). Note how individual knots forming 
the far-UV-bright ring show up in this picture. 



controversial: it could have formed together with the galaxy it- 
self, especially taking account that tidal stellar and gas arms are 
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often offset (e.g. Mihos 2001) . On the other side, though the 
projected image rules out that NGC 4262 is a classical polar ring 
galaxy ( i.e. the ring is not at 90 degrees to the galaxy) and the 
knotting of the ring is similar to the knotting instability seen in 
the not-polar-ring Cartwheel galaxy, the presence of a ring sur- 
rounding a SBO reminds shell and polar ring early-type galaxies, 
commonly thought to have undergone accretion and/or merging 
phenomena (e.g. Marino et al. 2009). 

To our goal we want (i) to show that both star-forming and 
HI rings surrounding this SBO galaxy (see Fig. 2) share the same 
radial distance from the galaxy center and spatial orientation and 
(ii) to model the kinematics of the ring(s) and of the galaxy body 
taking advantage of the existing innermost ionized gas velocity 
field (Sarzi et al. 2006) and the outer neutral hydrogen rotation 
curve (Krumm et al. 1985). 

2. Data and data reduction 

We de-archived both GALEX FUV (1344-1786 A) and NUV 
(1771-2831 A) images of NGC 4262 from deep frames (expo- 
sure time = 9.1 ksec [NUV], 7.2 ksec, [FUV], respectively) be- 
longing to the program GI2 017 (P.I. Zhong Wang). The tele- 
scope has a very wide field of view (1.25 degrees diameter) 
and a spatial resolution ~4".2 and ~5".3 FWHM in FUV and 
NUV respectively, sampled with 1.5"xl.5" pixels (Morrissey et 
al. 2007). 

We used FUV and NUV background- subtracted intensity 
images from the GALEX pipeline to compute the ultraviolet in- 
tegrated photometry and light profile of the galaxy. Background 
counts were estimated from the sky background image and 
high resolution relative response map provided by the GALEX 
pipeline itself. In addition we used optical Sloan Digital Sky sur- 
vey (SDSS) archival images (Adekman-McCarthy et al. 2008) in 
the g (3630-58830 A) and i (6340-8630 A) wavebands to obtain 
optical surface brightness profiles and (g-i) color profiles. 

Both UV and optical surface photometry have been derived 
by means of the IRAffl STSDAS ELLIPSE routine. ELLIPSE 
computes a Fourier expansion for each successive isophote 
(Jedrzejewski 1987), producing radial luminosity and geometri- 
cal profiles, sampled in nested ellipses, 2 arcsecs wide. These 
profiles are plotted versus galactocentric distance along the 
semi-major axis in Figs. 3 and 4. 

In order to estimate the errors on UV (AB) magnitudes, we 
propagated the Poisson statistical errors on source and back- 
ground counts. In addition to the statistical error, we added an 
uncertainty to account for systematic inaccuracies in the zero- 
point of the absolute calibration of 0.05 and 0.03 mag for FUV 
and NUV respectively (Morrissey et al. 2007). Our measured 
NUV total magnitude of NGC 4262 is consistent within errors 
with that measured by Gil de Paz et al. (2007). We derived also 
the FUV and NUV integrated fluxes of each individual UV- 
bright knot in order to estimate its current star formation rate 
(SFR). They were computed as m(AB) uv = -2.5 x log CRuv 
+ ZP where CR is the dead-time corrected, flat fielded count 
rate, and the zero points ZP= 18.82 and ZP=20.08 in FUV and 
NUV respectively (Morrissey et al. 2007). We corrected both UV 
fluxes and luminosities by adopting the foreground reddening 
given by Schleghel et al. (1998), namely E(B-V)=0.035. With 

1 IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO). NOAO is 
operated by the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the National Science 
Foundation. 



reference to the extinction curve of Savage & Mathis (1979), this 
translates into the extinction values A FUV =0.28 and A NUV =032 
respectively having adopted /l e //= 15 16 A and 2267 A for the 
two above filters. 



3. Galaxy and UV-bright clumps photometry 

The UV and optical surface brightness and color profiles of the 
galaxy body are shown in Fig. 3. The luminosity profiles of 
NGC 4262 at various wavebands indicate the predominance of 
a peaked distribution, that appears linear when plotted in r 1/4 , 
inside ~16 arcsec and of a exponential distribution outside this 
radius. So, data was fitted with a double-component curve made 
by a sum of a r 1/4 law and an exponential law, spanning all the 
range of observed luminosity profiles. A preliminary fit was per- 
formed with a single component only, inside and outside 16 arc- 
sec, to define the starting parameters, that are the effective radius 
r e and the apparent surface brightness at the same distance ii e . 
Then a model with both components was run, by adjusting the 
r e and fi e parameters for both luminosity laws with the goal of 
minimize the residuals between observed and calculated surface 
brightness. To perform the fit the very inner region of the pro- 
files, affected by the seeing (PSF= 2" in the optical and PSF= 5" 
in the UV), were excluded. The best fitting parameters are listed 
in Table 1 . 

Then, the total magnitudes of both components, bulge and 
disc were calculated and the resulting total galaxy magnitude 
for every waveband was obtained. These values also are listed in 
Table 1. We found a very good agreement with literature values, 
within a few hundredths of magnitude in i, g and NUV values. 
The error in effective radii and surface brightness were calcu- 
lated from the errors in the fitting parameters (slope and inter- 
cept) and are also reported in the table. Errors in total magnitudes 
were obtained from the correlation coefficient of the model fit vs 
the observed values. 

As one can see the disk contribution decreases as we go to- 
wards shorter wavelengths, and in the FUV profile is only barely 
visible, In the Fig. 3 an arrow marks the radial position where 
the disk starts to prevail. We note that — unlike optical color g- 
i — the FUV-NUV color becomes much bluer in the innermost 
region of the galaxy (the innermost UV color point represents 
the mean value inside the GALEX PSF). 

In Fig. 4 we show also the ellipticity and position angle pro- 
files, for all the four bands. The galaxy is seen almost face on, 
as the apparent disk ellipticity is 0.1 corresponding to an inclina- 
tion of ~26°, assuming an intrinsic flattening of 0.25 for the disk. 
In the g and i images the bar signature is clearly visible from 5 
to 20 arcsecs, with an increase of ellipticity (till e^ r =0.4) and a 
sharp change in the major axis to P.A=25°. This is indeed the po- 
sition angle of the major axis of the bar and is almost coincident 
with the P. A. of the major axis of the gas that reproduces its ob- 
served kinematics (both of the neutral and ionized components) 
as described in Section 5. In the NUV ellipticity and P. A. pro- 
files the signature of the bar is still present in the same region, 
whereas almost disappears in the FUV profiles where only the 
bulge component is visible. We note that the bar in the NUV is 
rounder than in the optical (ebaruv =0.3). In Fig. 5 (upper panel) 
is shown a schematic view of the galaxy (the apparent flattening 
on the sky of the stellar disk (black circles), the bar (red ovals) 
and outer HI ring (green ovals), while in Fig. 5 (lower panel) 
the edge-on view of the system with the line of sight is reported. 
On the upper right of each figure are given the used geometrical 
parameters. 
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Fig. 2. Comparison of optical (B+HST/ACS F475W) image of NGC 4262 (upper left), with GALEX NUV and FUV images (upper 
and lower right, respectively, and the smoothed HI column density (lower left). Contours are 0.8, 1..6, 2.4, 3.2, and 4.8 x 10 20 atoms 
cm -2 . The field of view of each panel is ~5'.3x5'.3. The ring physical size turns out to be -8.5 kpc while the galaxy body — as seen 
in the NUV — has a radius of only ~4 kpc. 



In addition to the luminosity profiles we measured the inte- 
grated mAB magnitudes and colors of the UV-emitting clumps 
that are given in Table 2. The brightest clumps we measured 
were selected by optical inspection and all the magnitudes are 
obtained inside an aperture of 7 arcsec in diameter. At the HST- 
based, adopted distance of 14.6 Mpc (Jordan et al. 2005) our 
(GALEX) spatial resolution is -400 pc and we cannot estab- 
lish whether their physical size is in line with typical Galactic 
loose OB associations (-80-100 pc). In fact the clump measured 
FWHMs are comparable with the large GALEX PSF (4-6") and, 
as such, they are not resolved. As stated above, the adopted dis- 
tance to convert AB magnitudes to luminosities is 14.6 Mpc. A 
map of individual measured, UV-bright clumps along the ring is 
shown in Fig. 6 Such a figure shows the full properly smoothed 
HI column density map (Krumm et al 1985) superposed to the 
recent GALEX image of the galaxy and its ring. The resulting 
match between the UV and HI rings is outstanding and clearly 
represent the same feature. What is more, the regions of highest 
HI condensation correspond to the brightest complexes in the 
UV. 



4. Star formation rate 

The present-day star formation rate of each UV-bright blob 
along the ring can be derived — following Kennicut (1998) — 
using its UV continuum luminosity and the relation 

SFIW (Mo yr _1 ) = 1.4 x 10- 28 L TOy (ergs s" 1 Hz" 1 ). 

FUV fluxes, luminosities and SFRs are given in Table 2. Note 
that our UV evidence of star formation within the galaxy body 
is only apparently in contrast with the lack of star formation de- 
rived in the infrared by Shapiro et al. (2010). In fact their low 
ratio between 8.0 jim and 3.6 jim SPITZER fluxes could be as- 
cribed to a low amount of dust as well. 



5. Galaxy and ring kinematics 

The innermost kinematics of NGC 4262 appears quite complex. 
Sarzi et al. (2006) write: "... this strongly barred galaxy shows 
an integral- sign pattern in the gas distribution and a twisted gas 
velocity field. In addition this object shows peculiar asymmet- 
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Table 1. Photometric parameters from the fit of the luminosity profile 
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5.2+0.92 
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NUV 
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Fig. 6. GALEX FUV image of NGC 4262 with superimposed the labels of individual, photometered UV-bright clumps along the 
outer ring. Their photometric parameters are listed in Tab. l.The field of view is ~4'.3x4'.3. The red empty circle shows the size of 
the photometric aperture used for all UV blobs (7 arcsec diameter), while the white empty square on the galaxy nucleus represents 
the size of the SAURON field of view (33" x 41"). 



ric distributions for the central values of sigma and for OIII/H/3 
ratio. The gas and stellar kinematics are strongly decoupled." 

To reproduce the gas shape and kinematics we used a model 
composed by a set of rings (described in Appendix A) inclined 
with respect to the sky plane and to the galaxy's reference plane. 
The intrinsic gas rotation curve is a typical Brandt (1960) curve 
with a steep increase in the first arcsecs from nucleus and an al- 
most flat part outside (Fig. 7). The observed values for atomic 
gas distribution and kinematics are from HI observations by 
Krumm et al. (1985) while the ionised gas kinematics has been 
studied by Sarzi et al. (2006). We deduced the values from the 
bidimensional maps and graphs of Krumm et al. (1985) and from 



a table of data used for the Sarzi et al. (2006) work and kindly 
sent us by dr. Marc Sarzi (the size of the SAURON velocity field 
is overplotted as a white rectangle in Fig. 6). For optical emis- 
sion lines error bars have been included as the r.m.s. of the ve- 
locities of each line. For HI points error bars have been taken 
from Krumm et al. (1985). The rotation curve has been extracted 
along a RA. of 29°, that corresponds to the sky elongation of the 
outer HI and UV ring, from 80" to 120". Then, a 100-ring model 
has been generated that fits the ring shape in the inner part (flat- 
tening and RA.). This generates the i, (p geometrical parameters 
described in Appendix A. The maximum extent of the model 
(100-th ring) has been adjusted to 120", the maximum HI ring 
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Fig. 3. UV/Optical surface brightness and color profiles of NGC 
4262 main body. In each panel we always plot the sum of the two 
single fitted profiles (arrows mark the radial distance where the 
disk starts to prevail). Upper panel shows the optical (g,i) pro- 
files, while mid panel gives the UV (NUV, FUV) corresponding 
profiles. As one can see the disk contribution decrease as we 
go towards shorter wavelengths, and in the FUV profile is only 
barely visible. Both optical (g-i) and UV (FUV-NUV) color pro- 
files are presented in the bottom panel. Note that the innermost 
UV color point represents the mean value inside the GALEX 
PSF. 



extent. As a first approach, no warping or twisting of the set of 
rings has been applied (in the model A£„=0 and Aa n =0 for all 
the rings). We found, in agreement with Krumm et al. (1985), 
that the inclination of the set of rings is 65° and their line of the 
nodes is atP.A.=29°. 

After a solution for geometrical parameter has been found, 
we tried different solutions to fit the observed rotation curve both 
of HI and ionised gas. As a result, we found that the maximum 
rotation of the ring holds at ~13" with a maximum velocity of 
~210 km/s ±10 km/s (see Fig. 7). The inner kinematics of the 
ionised gas may be influenced by the bar, that appears elongated 
at the same P. A. of the ring. The peaks of velocity that appear 
in [O III] and UJ3 at few arcsec fron nucleus may be due to this 
deviation from circularity. Unfortunately, no observed data for 
galaxy motions between 20" and 80" exist, so we tried the sim- 
plest possible solution fitting inner and outer data. 

The intensity distribution of [O III] and H/3 lines published 
by Sarzi et al. (2006) indicates a possible elongation of the gas 
to P.A.=0°. If real, this may indicate a twisting or spinning in 
the space of the gas ring toward the center. As an alternative, 
they may be a surposition of the gas circulating inside the bar in 
elliptic streamings with that of the disc in circular motion. But 
the map of the ionised gas distribution appears quite irregular, so 
we preferred do not search solutions with warped discs, to avoid 
increasing the number of variables present in the model. 



6. Origin of the stellar ring 

Taking into account the clear symmetry of the HI ring around 
NGC 4262, it is unlikely it comes from a primordial cloud of 
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Fig. 4. Position angle and ellipticity profiles for both optical and 
UV images of NGC 4262. Upper left panel: optical ellipticity 
profile. Upper right panel: P.A. optical profile. Note the remark- 
able change both in ellipticity and P.A. corresponding to the op- 
tical galaxy bar. Lower left panel: UV ellipticity profile. Lower 
right panel: P.A. UV profile. Note the disappearance of the bar's 
effect at ultraviolet wavelengths. 



neutral hydrogen such that discussed in detail by Thilker et al. 
(2009). On the contrary, Bekki et al. (2005) do model the strip- 
ping of rings and arcs of cold gas as due to interaction with other 
galaxies. In this context, a stream of gas pulled out of the disk 
of the galaxy (or contributed by a perturber) forms stars when 
it is compressed. Such an interaction scenario is supported also 
by a further hint, namely the observed inner decoupling of gas 
and stars velocity fields in NGC 4262 (Sarzi et al. 2006). In the 
specific case of NGC 4262 the existence of a mutual interac- 
tion with the other Virgo galaxy NGC 4254 has been proposed 
by Chyzy et al. (2007), while another example of such merger- 
induced, ringed lenticular galaxies could be NGC 404 (Thilker 
et al. 2010). As far as the UV-detected, disk's star formation is 
concerned, is likely occurring independently of the above possi- 
ble (ring-forming) interaction event. 



7. Summary and Conclusions 

Thanks to the UV- sensitive GALEX satellite we were able to 
detect an extended, UV-bright ring surrounding the otherwise 
normal SBO galaxy NGC 4262. Such a feature (not recogniz- 
able in the optical) appears to host several knots, likely con- 
sisting of hot star clusters. In this respect, NGC 4262 — having 
clustered UV-bright sources in its outer parts — could be clas- 
sified as a Type 1 extended ultraviolet disk (XUV). About the 
origin of such a structure, one should be aware that theoretical 
models (e.g. Bekki 2005; Higdon & Higdon 2010) ascribe the 
onset of rings and arcs of cold gas as well as the formation of 
young star rings to the past interaction with other galaxies. As a 
consequence, taking account also of the observed inner decou- 
pling of gas and stars velocity fields (Sarzi et al. 2006), we are 
pretty confident that a past major interaction episode underwent 
by NGC 4262 is responsible of the onset of the UV-bright ring 
we see today. 
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Fig. 5. A schematic view of the galaxy. Upper panel: the appar- 
ent flattening on the sky of the stellar disk (black circles), the 
bar (red ovals) and outer HI ring (green ovals) are shown. Lower 
panel: Edge-on view of the system with the line of sight reported. 
On upper right of the figures the used geometrical parameters are 
reported. 
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Appendix A: System geometry 

We define an observer reference system X"Y"Z" whose Y"Z" 
plane coincides with the sky plane and the X" axis represents 
the line-of-sight. The HI distribution is assumed to lie on the 
XY plane of a galaxy reference system in a set of concentric, cir- 
cular rings. The two sistems are inclined by the angles i and 
(see Figure lAH , that are the inclination and the azimuth of the 
line of sight with respect to XYZ. The relation between the co- 
ordinates of the two systems is described by the rotation matrix 
R, as discussed in Galletta (1983): 



(x"\ 

y 



'R\\ Rn R\3 
Rn R22 R23 
YR31 R32 R33J 



y 



(A.1) 



Each ring, labeled n, may be independent from the others 
and is inclined with respect to the galaxy reference plane by the 
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angles S n and a n , shown in Figure A.l. The change in S n repre- 
sent the warping of the gas plane, while the change in a n is the 
ring twisting, if present. The coordinate system of each ring is 
then 



7" 



(R" R" R 

lln Yin 
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\Zn) 



with the additonal conditions 



' = 
= 



Sky plane 
Ring plane 



(A.2) 



(A.3) 



The values of the coefficients R^ n (7, 0, S n ,a n ) are listed in 
Arnaboldi & Galletta (1993) and each ring is described by its 
radius r n and by an angle 0° < j3 < 360°. 



x n = r n cosp n 
y n = r n sinp n 



(A.4) 



To extract a simulated rotation curve, we calculate the cou- 
ples of value {y",z"} on the sky traced along a fixed RA. accord- 
ing to the relation: 



z = Ky 



(A.5) 



where K=tan(90+RA.). With this assumption, and the defini- 
tion of the ring plane (equation I A. 3 1) the equation [A]2] becomes: 



y f ; = r{R^ n cosp n +R^ n sinp n ) 

4' = Ky" = r(R^ ln cosp n + R^ 2n sinp n ) 



(A.6) 



that gives for each {y",z"} point a unique (5 value crossing the 
n-th ring: 



tgpn = ~ 



R 



3 hi 



K R'r!i 

2ln 



r>// v r>// 

K 32n ~ K K 22n 



(A.l) 



At this point, the rotational velocity of the n-th ring projected 
along the line of sight is 



Vossly", z\ n) = V circ (r n )(R'; 2n cos/3 n - R'{ ln sin/3 n ) (A.8) 

where \ rot is the intrinsic circular velocity of the rings. To 
interpolate the observed rotation values Y oss we used a simple 
rotation curve described by 



^circi^n) ~ 3 V m 



l+2(r n /r max ) 3/1 



(A.9) 



according to equation [26] and [28] discussed by Brandt 
(1960). 





X galaxy ^ 



Gnc eplan € 



slit\ 




Fig. A.l. top panel: The angles of the observer with respect to 
the galaxy reference plane, mid-panel: The inclination of a ring 
with respect to the galaxy reference plane, bottom panel: Image 
of the model rings plotted on a DSS I image of NGC 4262. A cir- 
cular "hole" has been created in the rings image only for show- 
ing purpose, to make the galaxy visible. The outermost ring is at 
200 arcsec, corresponding to the extension of the HI ring visi- 
ble in the other figures of this paper. The straight line across the 
rings indicate the position of the slit used for the extraction of 
the model velocity values. 
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Table 2. Integrated photometric properties 



Object 


R.A. 


DEC 


m FUV (AB) 


Ffc/v 


Lfuv 


m NUV (AB) 


Fnuv 




FUV-NUV 


SFR^f/y 




(ZVUU ) 


(ZUUU ) 


















Galaxy 


184.8776 


14.87777 


17.66+0.12 


2.466 


6.41 


16.26+0.01 


8.241 


21.43 


1.31 


8.9 


Clump 1 


184.8732 


14.89737 


22.68+0.39 


0.0308 


0.0800 


22.06+0.31 


0.0544 


0.142 


0.62 


0.112 


Clump 2 


184.8762 


14.88732 


21.18+0.24 


0.123 


0.318 


20.92+0.19 


0.156 


0.405 


0.26 


0.446 


Clump 3 


184.8603 


14.87558 


20.69+0.19 


0.192 


0.500 


20.59+0.18 


0.212 


0.548 


0.10 


0.700 


Clump 4 


184.8589 


14.86833 


21.61+0.31 


0.082 


0.214 


21.23+0.24 


0.117 


0.304 


0.38 


0.300 


Clump 5 


184.8554 


14.86131 


21.32+0.22 


0.108 


0.280 


21.15+0.23 


0.123 


0.327 


0.17 


0.392 


Clump 6 


184.8801 


14.86162 


22.13+0.31 


0.0510 


0.133 


21.70+0.30 


0.0759 


0.197 


0.43 


0.186 



Notes: (1) Coordinates are given in degrees 

(2) Fluxes are in units of 10~ 27 erg s _1 cm ~ 2 Hz~ l 

(3) Luminosities are in units of 10 26 erg s _1 Hz~ l 

(4) SFRs are in units of 10~ 2 M Q yr" 1 



